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TORSIONAL STIFFNESS OF THIN-WALLED SHELLS HAVING 
REINFORCING CORES AND RECTANGULAR, TRIANGULAR, 
OR DIAMOND CROSS SECTION 
By Harvey G. McComb, Jr. 


SUMMARY 


A theoretical investigation has been made of- the Saint-Venant torsion 
of cert ain composite bars. These bars are composed of two materials - 
one mater ial in the form of a thin-walled cylindrical shell and the other 
material in the form of a core which fills the interior of the shell and 
is bonded to it. 

An approximate boundary-value problem is formulated on assumptions 
similar to those of the theory of torsion of hollow thin-walled shells 
(Bredt theory). This boundary- value problem is solved exactly for a 
rectangular cross section and approximately for slender triangular and 
slender diamond cross sections. Results for the torsional stiffness 
constants are presented graphically. 


INTRODUCTION 


Certain airframe components such as wings, stabilizing fins, control 
surfaces, and helicopter rotor blades have been fabricated by employing 
a high-strength shell bonded to a core made of some light-weight material. 
The shell is formed in the external contour of the component and the core 
fills the interior and acts to stabilize the shell against local b uckling . 
Such a structure has been called a "foam-filled shell" because the core 
is often a foamed-plastic material. Metal honeycomb and balsa wood have 
also been used for cores. 

A large amount of literature exists on the problem of torsion of 
homogeneous isotropic cylindrical bars, but relatively little work has 
been done on the torsion of composite cylinders. A few exact solutions 
to problems in the torsion of composite sections are presented in refer- 
ences 1 and 2. Solutions for other cross-sectional configurations are 
desirable, and the methods used in references 1 and 2 do not .appear to 
,be applicable for sections of the type considered in the present report. * 
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In this report an approximate boundary-value problem for the tor- 
sion of foam-f ille d shells is formulated. The fact that the t hickne ss 
of the outer shell is small relative to overall dimensions of the cross 
section allo ws an approximation similar to that of the Bredt theory for 
the torsion of hollow thin- walled shells with free warping. (See ref. 3* 
pp. 298-302.) For a rectangular cross section, an exact solution to the 
approximate boundary- value problem is obtained. For slender tria n gul ar 
and slender diamond cross sections, approximate solutions which appear 
to be reas onab ly accurate are obtained. Finally, the results are com- 
pared with results based on an elementary concept of the torsional stiff- 
ness of foam-filled shells. 


SYMBOLS 


A,B 

arbitrary constants 

a,b,c 

cross-sectional dimensions (see fig. 4) 

C 

contour defining boundary of a region 

£U),g(£) 

arbitrary functions 

Gl,®2 

shear moduli of shell and core materials 

hn(x) 

Fourier coefficient (see eq. (21)) 

J 1>J2 

torsional stiffness constants 

K - Ggb - 

G 

m2 ~ 2 ~ 

Gqt 

Gjt 

kln' 1 ^ 

arbitrary constants 

M 

moment on cross section 

m 

index in equations (37) a^id. (6l) 

n 

summation index 


R 

s 

t 


region 

direction tangential to a contour in cross section 
thickness of shell wall. 
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to 

u 

u= GgU 


x,y,z 


7 

e = b/a 

e 

An 


dimension of cross section (see fig. 4) 
total complementary i energy per unit length 

coor dina tes along X-, Y-, and Z-axes, respectively 
arbitrary constants 
shear strain 

angle of twist per unit length 
eigenvalues 


\i = \/l + e 2 


V 

direction normal to a contour in cross section 


nondimensional coordinates in x- and y-dir e ct ions , 
respectively 

T 

shear stress 

<f 

stress function 

ft 

constant of integration 


Subscripts : 

1,2 refer to regions and contours indicated in figure 1 

Prime denotes differentiation with respect to the indicated 
variable. 


BASIC EQUATIONS 


In this section equations of elasticity are established for the 
Saint- Venant torsion of certain c ylindr ical bars composed of two materials. 
These equations are then particularized to the case of foam-filled thin- 
walled shells. Finally, the energy approach to the torsion of foam-filled 
shells is discussed briefly. 
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Torsion of Cylindrical Bars Consisting of Two Materials 

Consider a long cylindrical bar composed of two isotropic materials 
in which one material surrounds the other. A cross section of such a bar 
is shown in figure 1. For torsion with free warping, the stresses are 
given in terms of a stress function <f> as follows:' 


( 1 ) 


where i = 1 or 2. Each function ^ must satisfy Poisson's equation 

V 2 ^ = -2G±0 (2) 


in its corresponding region R^. 

The boundary conditions can be expressed in terms of the stress 
function by consideration of the components of shear stress normal and 
tangential to a contour. These components are, respectively. 


( 3 ) 


The condition that the external boundary of the section must be free of 
stress is obtained by integrating the first of equations (3) along C^: 


= ill = Constant (if) 

Cl 
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The conditions ■which must be satisfied at the interface C 2 can be seen 
by referring to figure 2. The tangential strain must be continuous across 
C 2 « In terms of the stress function', this condition is 


JL 5^1 

Gi dv 


>C 2 


_ 1 _ ^2 
G2 dv 



(5) 


Lastly, the shearing-stress component normal to C 2 must be continuous 
across C 2 , or 


When equation ( 6 ) 


c$2 

ds 


'C 2 


is integrated, the 



following equation results: 


( 6 ) 



(7) 


The problem is to find stress functions fa and fa which satisfy 

equation ( 2 ) in their respective domains and the boundary conditions 
(eqs. (4), ( 5 ), and ( 7 )). 

The total moment on the cross section is given by 

dRi (8) 



This equation can be written in terms of 0 as follows: 
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Integrating by parts 


and making use of equations (4) and (7) gives 


M = y ftjL F (-x dy + y dx) 
i=l J C± 



( 10 ) 


The stresses and moment on the cross section are inde p endent of 

therefore, these constants may be chosen arbitrarily. For convenience, 
ftl and ft2 are both taken to be zero so that the expression for the 

moment is analogous to that for the torsion of a uniform ctosb section; 
that is, moment equals twice the volume under the ^ diagram. 


Torsion of Foam-Filled Shells 

The equations of elasticity are particularized for the case of a 
cylinder made of a thin-walled shell of one material filled with a core 
of another, that is, a foam-filled shell. In figure 3 a general cross 
section for such a c ylinde r is illustrated. Because the thickness of 
the shell wall is small compared with the overall dimensions of the ctosb 
section, the stress in the wall can be assumed to be uniformly distributed 
over the thickness. This stress is equal, to the normal derivative and is 
given by 


‘wall “ 


d01 _ 01 

dv t 


Q2 


( 11 ) 


where T t wl .n represents the stress in the shell wall. With the use of 
equation (ll), equations (^), ( 5 ), and ( 7 ) can now he written as 


and 


01 



= 0 


d0 2 

dv 


_ _ °2 01 

G 2 02 

Gx t 

Gq t 

C 2 

C 2 


C2 


( 12 ) 


( 13 ) 
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Let the contour C in figure 3 he the middle surface of the shell 
wall. If the shell vail is assumed to be concentrated at its middle sur- 
face, then C can be thought of as representing the interface, middle 
surface, and outer boundary of the w all . Consider R as be ing the region 
bounded by C. The problem can now be formulated as follows. Find a 
function 0 satisfying the equation 

= -2G20 • (l4) 


in R and the equation 


dv 



= _G2 g 

t 


( 15 ) 


along C. The moment on the cross section is equal to twice the -volume 
under the <f diagram, or 


M = 



dx dy 


(16) 


Energy Approach 

Approximate solutions for the torsion of foam-filled shells can be 
obtained by the energy method. The complementary energy for such a body 
is the sum of the stress energy of the core, the stress energy of the 
shell wall, and the negative of the work done by the external moment 
acting through the angle of twist. For the composite cross section 
shown in figure 3 the complementary energy per unit length is 


U = i + T yz 2 ) ^ ay + -^rf c ( T wanft is - mg (17) 

R 
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where r xz and Ty Z represent the shear stresses in the core and. r ^n 

represents the shear stress in the shell wall. In terms of the stress 
function 0, U becomes 

n - 4^ (k 2 + ** ' ^ * + l *■ (18) 

R 

where the subscripts on 0 denote the partial derivative with, respect 
to the indicated variable. 

When the variation of U is equated to zero and integrations by- 
parts are carried out, the following equation is obtained: 



(19) 


It is seen that, if 60 is arbitrary in R and along C, equations (l4) 
and ( 15 ) must be valid. 


The torsion of a foam-filled shell is analogous to the problem of 
the deflection of a membrane stretched over the region R, subjected to 
lateral pressure, and supported along the contour C by infinitely many 
springs which are constrained to distort only in the direction normal to 
the plane of R. Some discussion of approximate solutions of problems 
of this type is given in reference 4. 


SOLUTION OP SPECIFIC PROBLEMS 


In this part of the report an exact solution for a rectangular cross 
section is obtained by satisfying equations (l4) and ( 15 ) . For the slender 
triangular and slender diamond cross sections, exact solutions do not 
appear feasible, and approximate solutions are obtained instead by using 
equation (18) . Two appro ximat e procedures are utilized in each case; 
the first is the Rayleigh-Ritz method and the second is a more general 
variational procedure, herein called the "variational method. " This 
latter procedure is applied, for the most part, in cases where the 
Rayleigh-Ritz method becomes cumbersome. These two approximate methods 
may be applicable to other sections of practical interest for which solu- 
tions are not available. 
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Rectangular Cross Section 

The notation for a rectangular cross section is shown in figure 4(a). 
The thickness of the shell wall is assumed to be constant. The problem 
is to f i nd a function 0 which satisfies equation (l4) within the rec- 
tangle and the following boundary conditions: 


_ _ g 2 £ 

(at x a + a) 

(20a) 

dx + Gi t 



M 

(at y = ± b) 

(20b) 

dy + G]_ t 



The function <f> is, of course, symmetric about both the X- 

and 

1 


7-axes. The symmetry condition about the X-axis is satisfied when 
is taken in the form 

00 

$ - 21 bn( x ) cos Any ( 21 ) 

n=0 

where the functions cos Tv^y form an orthogonal set. The right-hand 

side of equation (l4) can be expanded in a Fourier series of the func- 
tions cos Any in the interval -b ^ y ^ b, and this expansion yields 


V20 - - 


4 G 20 sin A^b 


n?=0 Anb + sin A^b cos A n b 


cos Any 


( 22 ) 


Substituting the assumed solution (eq. (21)) into equation (22) and 
equating coefficients of like terms gives the following ordinary dif- 
ferential equation for h n (x) : 


4G28 sin A^b 


hn(x) " - A n 2 hn.(x) = - 


Ajjb + sin Anh cos An.b 


(23) 
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The solution to equation (23) is 


4 G 2 0 sin Xjjb 

h n (x) = k-^ sihh A n x + k^ cosh A n x + — 

An (Anb + sin Anb cos Anb) 


(24) 


The constants k-|n vanish because of symmetry. The constants and. 

the eigenvalues An can be found from the boundary conditions at x = a 
and y = b, respectively. Consider first the condition at y = b. The 
substitution of equation (21) into equation (20b) yields 


Go 

An sin A n b = — - cos A n b 

Git 


(25) 


Therefore, the eigenvalues are given by 


tan Anb = 


K 

Anb 


(2 6) 


where 


Gob . 

K = 

G 1 t 


At x = a, the substitution of equation (21) into equation (20a) yields 


kgnV sillh V 


°2 

G]t 


k2n cosh Ana + 


4 G 29 sin Ajjb 

An 2 (Anb + sin A n b cos A n b) 


( 27 ) 
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Therefore, 


&2 n = 


-4EG 2 9b 2 sin T^b 


A n 2 b 2 (A a b slnh An.a + K cosh A n a) (Anb + sin Anb 008 ^n b ) 


( 28 ) 


Consequently, the stress function is 


0 = 


00 


te20b 2 sin Ajjb cos A^y 


1 - 


K cosh AjjX 


n=0 A a 2 b 2 (A n b + sin A n b cos A n b) \ A n b sihh A a a + K cosh A n a. 


(29) 


The moment on the cross section is given by the formula 


cbc dy (30) 


in terms of either or Gg 


1=0x1!= C&2 (3D 


The torsional stiffness constants and Jg are obtained upon sub- 

stituting equation (29) into equation (30) and carrying out the indicated 
integration. When the results axe expressed in the form of equation (31), 
it is seen that 

\ 

Jl = ij-ctQ^KA 


M 


- e/T „ 

J o J o 


The torsional stiffness can be expressed 
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where 


sin 2 Anb 


it=0 An^b^CAnb + sin A n b cos Anb) 


K sinh Ana 

Anb(Anb sinh Ana + K cosh Ana) 


The series A converges very rapidly. For various values of K, the 
eigenvalues are easily located from the intersections of the hyper- 
bola K/Anb and the curves tan A n b as indicated in figure 5 * Plots 
of Jx and Jg against the cross-section aspect ratio to/c for various 
values of K are presented in figure 6. 


As the stiffness of the shell-wall material vanishes, K 
infinity and A n b approaches - 2 * 1 - — 1 - — . It can be shown that the 


approaches 


limiting value of J2 is the solution obtained by the theory of elastic- 
ity for the torsion of a homogeneous rectangular ctobs section as given 
in reference 3 (p. 278). 

As the stiffness of the core material, vanishes, K approaches zero 
and, the solution should, approach that of the Bredt theory for the torsion 
of hollow thin-walled shells with free warping. As K approaches zero, 

Anb approaches me. It is obvious, then, that all the terms in A vanish 
for which n > 0 . Investigation of the term for which n =» 0 shows 
that Jx does approach the result given by the Bredt theory which is 


lim Jx = 
K — >0 


2ct 0 2 t 




where Aq is the area enclosed by the median line of the shell wall. 


Slender Triangular Cross Section 

If the energy approach is used, two approximate solutions are 
obtained for a cross section in the shape of a slender isosceles triangle *■ 

with a constant-thickness shell wall as illustrated in figure 4 (b) . One 
solution is obtained -by the Rayleigh-Kitz method and another solution is 
found by utilizing the calculus of variations and the boundary- layer •> 

technique of reference 5* 
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Rayleigh-Ritz method .- In terms of nondimensional coordinates , 
complementary energy (eq. ( 18 ) ) , can be -written for the triangle as 


the 

follows: 


U= G 2 U 




i dtj d£ + 




(35) 


where 


I = 


x 

a 



p = |1 + e 2 


Note that the equation of the sloping side of the triangle in nondimen- 
sional coordinates Is simply rj - £. 

The stress function 0 must be an even function through the thick- 
ness, and for slender sections it is usually sufficient to assume a par- 
abolic variation in the thickness direction. For thick sections, how- 
ever, it may be necessary to include additional terms in the thickness 
direction to get satisfactory accuracy. Suppose now that 0 is assumed 
to be a polynomial 


0 - 5Z On£ n + 0 ti 2 (36) 

n=0 

When equation ( 36 ) is substituted into equation (35) 8 E& U Is m i nim ized 
with respect to the parameters c% and 0 , the following r + 2 equa- 
tions result : 
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2k 



where 0 ^ m ^ r. 

Solution of the system of simultaneous equations (37) yields 
and p. From equation (l6), the moment, is 
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When K is large compared with unity a large number of equations may 
be required for reasonable accuracy. An approach which avoids this 
difficulty is developed in the succeeding section. • 

Variational method .- Instead of assuming for $ the polynomial of 
equation (36), suppose arbitrary functions of | are allowed to remain 


equation \5b) , suppose arbitrary an 
and 0 is taken to be of the form. 


f(S) + il 2 g(!) 


When_e quat ion (40) is substituted into equation (35) and the variation 
of U with respect to admissible variations in f and g is equated 
to zero, two simultaneous ordinary differential equations for the func- 
tions f and g are obtained as follows: 


e 2(£f*)« + ^i3g*)’ - Kji(f + |2g) = -2G20b 2 | 


£(|3f *) * + ^(i5g«) 1 - K^(.|2f + |4g) - | |3g - - | G20b 2 |3 
3 5 3 3 


and the following boundary conditions .are obtained: 


■M' 1 ) 


(42a) 


^|3 f « + | |5 g 


(42b) 


(tf ' + Kf) = 0 


(42c) 
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eg' + Kg ) =0 (42d) 

k '£=1 


where the primes denote derivatives with respect to g. 

The differential equations (4-1) are linear with variable coeffi- 
cients, and. it appears to be a difficult task to find, an exact solution 
to the system. For the case of slender cross sections, however, an approx- 
imate solution is possible by utilizing the "boundary layer" technique 
discussed in reference 5» Notice that the differentiated terms in equa- 
tions (kl) are multiplied by e 2 , a quantity which for slender cross 
sections is small compared with unity. Differential equations having the 
most highly differentiated terms multiplied by a small quantity are char- 
acteristic of the type of b oundary- layer problems considered in 
reference 5 . 

Suppose, initially, that f and g are slowly varying functions 
throughout the region 0 ^ | 1. The term "slowly varying 1 is intended 

to mean that the maximum values of the functions f and g and their 
derivatives which appear in equations (kl) are of the same order of magni- 
tude. Then, as long as K is at least of the order of unity, the terms 
in equations (kl) which contain e 2 have little influence on the solu- 
tion. Consequently, a good approximation to a particular solution is 
obtained by Ignoring the terms in equations (kl) which contain e 2 . 

When this procedure is carried out, it is found that the approximate 
particular solution satisfies the boundary co n d it ions at g = 0 but does 
not satisfy the boundary conditions at g = 1. It can be concluded that 
the required solution is such that the functions f a n d g are not 
slowly varying everywhere in the region 0 < g 1. Somewhere the deriv- 
atives which appear in equations (kl) must take on values which are of 
the order of e~ 2 so that the terms containing e 2 can have an appre- 
ciable influence on the solution. 

It is assumed that the region where the derivatives of f and g 
have values of the order of e -2 is confined to a so-called boundary 
layer in the neighborhood of g = 1. On the basis of this assumption, 
the particular solution alone is a good approximation to the exact solu- 
tion away from g = 1. Then, by focusing attention on the boundary layer 
close to g = 1, it is possible to obtain an approximate homogeneous 
solution to equations ( 4l) which modifies the particular solution in 
such a manner that the boundary conditions at g = 1 can be satisfied. 

It is convenient to get a particular solution as a power series 
in e instead of ignoring completely the e 2 terms in equations (kl) . 
Assume that a solution can be expressed in the form 


v 
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fp(l) = fpo(l) + efpi( 6 ) + e 2 fp 2 (i) + • • • 
gp(i) = gpo(S) + £gpi(i) + s^gppCi) + • • • 


l (^3) 


where the subscript P denotes a particular solution. When equations (43) 
are substituted into equations (4l) and. coefficients of like powers of e 
are equated, pairs of simultaneous equations are obtained for the coef- 
ficients in the power series. For example, when the coefficients of the 
zeroeth power of e are equated there results 


K[i(fpo + S 2 gpo) = 2G20b 2 g 

► 

Kix(s 2 fpo + £^gpo) + I'^gpo = J G20b 2 |5 


(44a) 


When coefficients of the first power of e are equated the result is 


f Pl + § 2 gpi = 0 

Kp^^pi + i^gpij + ^ i^Spi = 0 


j 


(44b) 


When coefficients of the second power of e arfe equated the follow in g 
equations result : 


(l f P0') + I^Spo') ~ Kn( f P2 + £ 2 Sp2) = 0 
j^Po') + 5^ 5g Po' ) ~ K + l (^ 2f P2 + ^ g P2 ) " 3 &5g P2 


>■ (44c) 


0 


Si m ilar equations are obtained when coefficients of higher powers of e 
are equated. 
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Solution of equations (44) results in the following expressions 
for fp and gp: 


fp = G20b £ 


St + |2 + e 2/_2_ + il + S 2') + 

^2 ■ 141 


gp = -G20b 2 (l + e 2 + . . .) 



(45) 


Coefficients of the odd powers of e vanish. This particular solution 
satisfies the boundary conditions at g = 0 but not at | = 1. 

A homogeneous solution can be obtained which modifies the particular 
solution in the vicinity of | = 1 in such a way that the boundary con- 
ditions at g = 1 can be satisfied. In order to determine the homo- 
geneous solution it is convenient to put equations (4l) into a form in 
which, in the neighborhood of | = 1, the terms containing derivatives 
are of the same order of magnitude as the remaining- terms . Such a con- 
version is provided by the coordinate transformation 


g = 1 + eg 


C«!) 


When the transformation (eq. (46)) is introduced into equations (4l) and 
the right-hand sides are set equal to zero, the following equations are 
obtained: 


|(l + ef)f ’J * + |j(l + eI)Vj* - + (1 + el) 2 gj - 0 

i[(l + ef)V~| + if(l + ef) 5 g’l - Kjtfil + ef)^ + (l + el) k g 


- >( 47 ) 


|(1 + eg) 3 g = 0 

5 


where the primes now denote differentiation with respect to g. 
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The homogeneous solution can be expanded in powers of e 


%(!) = %)(*) + e^Hld) + e 2 fH2(S) + • • • 
gn(i) = Ssoii) + eg Tn ( I ) + g 2 gg-p(i) + • • • 


m 


» 


where the subscript H denotes a homogeneous solution. If these expres- 
sions are substituted into equations (lf-7) and the coefficients of like 
powers of e are equated, systems of ordinary differential equations 
are obtained. For the terms fjjg and gjjQ,. a set of homogeneous equa- 
tions is obtained: 


f Ho" + ^ sed" ~ + Sho) = 0 

j f Ho" + j ®Ed" “ + s Hd) ' j g H0 = 0 


m 


For the terms f m and gg^ the following nonhomogeneous equations are 
obtained: 


f Hi" + ^ skl" ■ + skl^ = ■^hd*) “ ^sho * ) + 2k + 1 £sho 

\ f m." + 5 ®hi" ~ + ~ \ s hi = ~(^ f Ho') " (^ g Ho’) + * 

2Kni ^f H0 + -2g2o) + ^ g HD 


Additional sets of no nho mogeneous equations would result for the coef- 
ficients of higher order terms. 
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It is found that neglecting terms of the order of e in equa- 
tions (48) is equivalent to neglecting terms only of the order of 
in the final result for the torsional stiffness constants. Therefore, 
a final result which includes an terms linear in e can be obt ain ed by 
solving only equations (49) and dropping all higher order terms in the 
homogeneous solution. Solutions to equations (49) are of the form 

ffiD = 



Substitution of these functions into equations (49) yields 



This system of equations has a solution only if the determinant of the 
coefficients vanishes. When the determinant is set equal to zero a 
biquadratic equation for A is obtained. The solution to the biquadratic 
equation is 


A = ± 


15 + 6*41 ± \] 


(15 + 6 Kji) - 60Kp. 


(52) 


The homogeneous solution must vanish as the absolute magnitude of | 
increases, and since f is always negative or zero, only the two positive 
roots for A are required. 

From the first of equations (51) the B coefficients can be written 
in terms of the A coefficients 


B = -3 



A 


A 2 - 3Kp 


' (55) 
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When the | coordinate is transformed to the £ coordinate anri terms of 
the order of e 2 are dropped, an approximate general solution to equa- 
tions (4l) is obtained 


f = ZI *ne e e"^ + G20b 2 ^ + S 2 J 

► 

2 > 2 _ Tfii An. 

g = ~ y~ 3 — %e“ 6 - C^Sb 2 

ft " 1 A n 2 - 3K+i 


(5^) 


This solution satisfies the boundary conditions at £ = 0. The boundary 
conditions at £ = 1 are used to determine Al and 

When equations (5*0 are used, the stress function <f> becomes 


4 = G20b 2 


An. 1 - 3 


rt=l 



-_***_ «2 


(55) 


A n - 3Kp 


where 



An 

G20b 2 


and where terms of the order of e 2 are dropped. The moment is cal- 
culated by substituting 4 into equation (l6) and performing the integra- 
tion. The results for the torsional stiffness constants and J 2 


are 
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where 



j 


r = 



6 




(56) 


(57) 


and where terms of the order of have been dropped. It is seen now 

that In calculating the arbitrary constants A n , terms of the order of e 
may be neglected. When equations (5*0 are substituted into the boundary 
conditions (eqs. (42c) and (42d)) the following expressions for A n are 

obtained after dropping terms of the order of e: 


Al» 


(l + K ) 


A 2 2 - 2K 


a 2 = - 


(Ai + K)^(15 + 6K+i) 2 - 60Kp 

(; ts k 2 - 5 


(58) 


(7)2 + k)^(15 + 6Kii) 2 - 60Kjj, 


Tn the limiting case where K approaches zero the boundary-layer 
technique becomes invalid. Therefore, the solution cannot be expected 
to approach the proper result for a hollow thin-walled shell. When K 
approaches infinity, an approximate solution for a solid cross section is 
obtained; and r is given by 


lim T = 
K — 


1 



e 


(59) 
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Slender Diamond Cross Section 

Rayle igh-Rit z method .- For the slender diamond cross section the 
notation is shown in figure 4(c). The complementary energy becomes 


U = 




i dT] ag + 2 

e 



< 1 ! 


(60) 


Substitution of the polynomial (eq. (36)) into equation (60) and minimiza- 
tion with respect to <% and p yields the. r + 2 equations 



where 0 <? m g. r. 


The moment and the torsional stiffness constants and J 2 are 

found by utilizing equation (l 6) . The expressions for the torsional 
stiffness constants turn out to be precisely the same as for the tri- 
angular cross section given in equations (39)* Of course, for the diamond 
cross section, c% and p are obtained from equations (6l) . 


Variational method . - Through the use of the calculus of variations 
and expression (40) for j , the differential equations (4-1) are found 
to be valid also for the diamond cross section. The boundary equations 
at | = 0 (eqs. (42a) and (42b)) also hold, for the diamond cross section. 
The boundary conditions at i = 1, however, are now given by 


f’ 


g 


= 0 


6=1 


= 0 


6=1 


(62) 


where the primes denote differentiation with respect to |. 
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The bou n d ar y-layer technique yields a general solution of the same 
form as that obtained for the triangle. However, now it is found that 
neglecting terms of the order of e in the homogeneous solution 
(eqs. (48)) is equivalent to neglecting terms only of the order of 
in the final expressions for the torsional stiffness constants. Conse- 
quently, for the diamond cross section a result which includes an terms 
of the order of can be obtained by solving equations (49) and dropping 

all higher order terms in the homogeneous solution. It is consistent, now, 
to keep the e 2 term in the particular solution. 

The approximate general solution to equations (4l) for the diamond 
cross section becomes 



The arbitrary constants Ax and Ag are determined by substituting 
equations (63) into equations (62) . 

The stress function becomes 


0 = G^b 2 




, v 2 - Kfi 

0 — 5 

A n 2 - 3Kp 



\i ^n. 

e" + 


2| 

Kji 


+ 


i 2 


T ] 2 + 





(64) 


If this expression for (f is substituted into equation (l6) and inte- 
grated, the moment can be calculated. The torsional stiffness constants 
become 



NACA TN 5749 


25 


tM 


ctQ^b 

Jl = 2KT* 

12 


et 

j 2 = J2L r* 
12 


where 


(65) 


r*= 1 + JL+ e 2(_6_ + 1.4. a. ^ 7 e / -2Ku 

* 


(66) 


When equations (63) are substituted into the boundary conditions 
(eqs. (62)) and terns of the order of e£ are neglected, the following 
expressions are obtained for the arbitrary constants: 


Ai = e 




2Xg + 3Kp 


^ ^1^15 + 6K41) 2 - 60Kn 


A 2 = - e 1 * ^1 + 3Kp 


% A 2 ^(15 + 6Kh) 2 - 6OK4J. 


(67) 


In the cane where 
for a solid section is 


K approaches infinity an approximate solution 
obt ain ed, and r* is given by 


lim T* = 1 _ 3 e 2 

K — >00 


( 68 ) 
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RESULTS AND DISCUSSION 
Presentation, of Results 


The results of the calculations outlined in this report are presented 
in figures 6, 7, and 8. These figures show plots of the torsional stiffness 
constants against cross-section aspect ratio for various values of the 
parameter K. The torsional stiffness can be expressed either in terms 
of the shear modulus of the core material or the shear modulus of the 
shell-wall material. The torsional stiffness constant associated with 
the shear modulus of the core materia]. j£ is plotted in parts (a) of 

figures 6 , 7, and 8, whereas the constant associated with the shear 
modulus of the shell-wall material is plotted in parts (b) of these 

figures. 

In figure 6 are shown the results of the exact solution of the dif- 
ferential equation (l4) with the boundary conditions (20) for a rectangular 
cross section. In figures 7 and 8 are shown the results of the approximate 
solutions (the Rayleigh-Ritz method and the variational method in conjunc- 
tion with the boundary- layer technique) for the triangular and diamond 
cross sections. A five parameter polynoMal was -used in the Rayleigh- 
Ritz method. 


Accuracy of Approximate- Methods 

Solutions by the Rayleigh-Ritz and variational methods also were 
obtained for the rectangular cross section. A comparison of these results 
with the exact solution provides an indication of the accuracy of the 
approximate methods. 

A polynomial with only three parameters was used for the Rayleigh- 
Ritz method in this comparison. The results showed that for all aspect 
ratios and for K less than about unity the stiffness given by Rayleigh- 
Ritz method is less than 3 percent lower than the exact stiffness. It 
is believed that the five-parameter Rayleigh-Ritz method used for the 
triangular and diamond cross sections should yield slightly more accurate 
results. Of course the accuracy of the Rayleigh-Ritz method can be 
improved for large values of K by including more terms in the polynomial 
for 0. However, the number of simultaneous equations which must be solved 
increases with the number of unknown parameters. 

The stiffness calculated by the variational method for the rectangular 
cross section was less than 1 percent in error for values of K greater 
than about unity and the aspect ratio to/ c less than about l/4. The 
boundary-layer technique yielded a slightly more accurate solution to the 
differential equations obtained for the rectangular cross section than for 
equations (4-1) which arise for the triangular and diamond cross sections. 



NACA TN 37^9 


27 


Thus the results of the variational method for the triangular and diamond 
cross sections are probably not quite as accurate as for the rectangular 
cross section. It appears that for slender cross sections (small values 
of to / c) the variational method is more accurate for large values of K 
and the Rayleigh-Ritz method is more accurate for small values of K. 

In figures 7 and 8 the solid curves are results of the five-parameter 
Rayleigh-Ritz method, and the dashed curves show results of the variational 
method. The Rayleigh-Ritz method certainly leads to a lower bo und for the 
torsional stiffness. The variational approach also leads to a lower bo und 
provided the boundary-layer technique yields a sufficiently accurate solu- 
tion to equations (^-l) . From the results on the rectangular cross section 
it appears reasonable to assume that for K ^ 1 the variational method 
probably gives a lower- bound. Therefore, for any particular case wher e 
K £3 1, the higher of the two values of torsional stiffness calculated by 
the two approximate methods is the more accurate. 


Remarks on an Elementary Calc ulat ion 


It is conceivable that a first approximation to the stiffness of a 
foam-filled shell could be made simply by adding the torsional stiffness 
of the core alone to that of the shell alone and neglecting the stiffening 
effect which results from bonding the core anfl shell together. Calcula- 
tions were made by this elementary procedure anfl the results are shown 


«1 eiem 

in figure 9* The ratio — is the ratio of as calculated by 

J l,th 

the elementary procedure (neglecting the bonding effect) to as cal- 


culated by the theory of this paper. It is seen that for the rectangular 
cross section the result of the elementary calculation is never more than 
5 percent lower t han the exact solution. For the triangular 
cross sections, however, the elementary procedure does not lead to such 
good results, and the discrepancy can be as much as 25 percent. 


CONCLUDING REMARKS 


An approximate boundary- value problem is set up for the Saint-Venant 
torsion of cylindrical thin-walled shells bonded to a core which fills 
the interior of the shell and which is made of a material different from 
that of the shell wall. Solutions for the torsional stiffness are obt aine d 
for three particular cross-sectional shapes - an exact solution to the 
bou n da r y- value problem for rectangular cross sections and approximate 
solutions for slender triangular and slender diamond cross sections. 

The approximate solutions are obtained by the use of two energy procedures. 
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These methods may be applicable to other cross sections of practical 
interest. The choice of approximate method for any particular problem 
depends on the range of parameters involved. 

The stiffness obtained by the simple procedure of adding together 
the individual stiffnesses of the core and the hollow shell (neglecting 
the effect of the bond) yields results less than five percent low for 
rectangular cross sections. For slender triangular and diamond cross 
sections this elementary approximation is generally not so good and in 
certain cases it yields results which are considerably low. 


Langley Aeronautical Laboratory, 

national Advisory C ommit tee for Aeronautics, 
Langley Field, Va., June 7, 195^. 
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Figure 1.- Cross section of cylindrical bar composed of two materials . 



Figure 2.- Shearing stresses on an element at the interface between the 
materials. Superscripts on symbols correspond to regions indicated 
in figure 1. 



Figure 3 .- Cross section of a composite thin -walled cylindrical shell 
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(b) Triangular cross section . 
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Figure k.- Notation used in analysis. 
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Determination, of eigenvalues in exact solution for rectangular 

cross section. 
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(b) Constant associated with shear modulus of shell -wall material 


Figure 6.- Torsional stiffness constants for a composite thin-walled 
cylindrical shell of rectangular cross section. M/9 = GiJi = GqJq/ 

n, tO 
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(a) Constant associated with shear modulus of core material . 



Figure Torsional stiffness constants for a composite thin- 
walled cylindrical shell of slender triangular cross section. 


M/0 = GiJ x = G 2 J 2 ; K = 
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